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We present a method for inverting charged particle velocity map images which incoorporates a non-uniform 
detection function. This method is applied to the specific case of extracting molecular axis alignment from 
Coulomb explosion imaging probes in which the probe itself has a dependence on molecular orientation 
which often removes cylindrical symmetry from the experiment and prevents the use of standard inversion 
techniques for the recovery of the molecular axis distribution. By incorporating the known detection function, 
it is possible to remove the angular bias of the Coulomb explosion probe process and invert the image to 
allow quantitative measurement of the degree of molecular axis alignment. 


I. INTRODUCTION 

Photofragment imaging has become a standard tool in 
the chemical physicists arsenal of tools for detailed mea¬ 
surements of processes in gas phase molecules 1 21 since its 
inception in the 1980^ and the subsequent evolution of 
the technique known as velocity map imaging (VMI)P 

In a photofragment imaging experiment an expand¬ 
ing sphere of charged particles is projected onto a two- 
dimensional (2D) position sensitive detector. The aim 
of this detection method is to extract the original 3D 
distribution of the charged particles from the 2D pro¬ 
jection, and so recover the energy and angular infor¬ 
mation regarding the fragmentation process. While it 
is possible to arrange to image only a central 2D slice 
of the 3D distribution using laser or electrostatic slicing 
technique^ ® it is frequently the case that experimental 
considerations require ima ging o f the 2D projection to 
regain the 3D distribution. 1 7: “ -21 Reconstructing the 3D 
distribution from a single 2D projection requires that the 
experiment is cylindrically symmetric about an axis ly¬ 
ing in the plane of the detector; if this condition is not 
met, then in general it is only possible to reconstruct 
the 3D distribution tomographically from multiple 2D 
pro j ections E3H1 

However, there exists another class of experiments in 
which the experiment lacks cylindrical symmetry due to 
a non-uniform detection function that is otherwise mea¬ 
surable or known. In these circumstances we show here 
that it is possible to invert a single 2D image to recover 
the 3D distribution while simultaneously correcting for 
the non-uniform detection function. 

Our methodology is motivated by the desire to char¬ 
acterize the degree of molecular axis alignment and ori¬ 
entation produced through the interaction of a molec¬ 
ular sample with intense non-resonant laser fields E3EH 
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For such aligned/oriented samples to be useful for 
subsequent experiments, such as x-raj® 21 or elect ron 
diffractionpSHMi high-order harmonic generation I 25 * 26 ! and 
photodissociation or photoionization studiespBiia it j s 
necessary to characterize and quantify the degree of 
alignment and orientation produced. It is common prac¬ 
tice to utilize Coulomb explosion imaging (CEI) for this 
characterization 13 15 11 19 where an intense probe laser 
pulse with duration much shorter than molecular rotation 
is used to remove multiple electrons from the molecules 
under study. The multiply ionized molecules subse¬ 
quently undergo rapid fragmentation due to Coulomb re¬ 
pulsion, and imaging the resulting ion fragments is then 
used to establish the orientation of the molecules prior 
to ionization; under the assumptions that the fragmenta¬ 
tion happens rapidly with respect to rotation, and that 
the fragments recoil in the direction of molecular bonds, 
there is a direct correlation between the fragment recoil 
and the molecular orientation. While this technique has 
been very succesful at analysing alignment and orienta¬ 
tion of molecular samples it has proved difficult to fully 
quantify the degree of alignment/orientation since the 
ionization process in CEI has a strong dependence on the 
molecular orientation with respect to the ionizing laser 
polarization. In fact it has become common practice to 
characterize the molecular axis alignment and orientation 
in such experiments using expectation values calculated 
for the resulting VMI image, such as (cos 2 $2d), where 
$ 2 d is the angle in the plane of the detector measured 
from the axis of laser polarization. This expectation 
value includes anisotropy due to the CEI probing, and 
is calculated in lieu of a method suitable for extracting 
the true molecular axis distribution. Here we show that 
it is possible, under certain circumstances, to remove the 
effect of the orientation dependence of the CEI probe 
from the measured molecular axis distribution, and so 
extract the true moments of the axis distribution. This 
is possible by making an independent measurement of the 
CEI orientation dependence using an isotropic gas under 
the same conditions as the alignment measurement which 
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FIG. 1. Relationship between coordinates of the original 3D 
Newton sphere (right) and the 2D projected image (left). 


can subsequently be deconvoluted from the CEI images 
of aligned/oriented samples. 


the position of a charged particle on the Newton sphere 
which is projected onto the detector (see Fig. [l]). We 
show in Appendix [A] that the Abel transformation may 
be written in spherical polar coordinates as 


F(R,e) 



Vr 2 - R 2 

rg{r,e)D(r,9,4>) ^ 
Vr 2 - R 2 r ' 


( 2 ) 


where F(R,Q) is the projected (image) distribution, 0 
is the polar angle measured in the detection plane with 
respect to the z axis, and R is the distance from the 
image centre (see Fig. [l]). 

Since Eq. (2| has a similar form to Eq. Jl]) many of 
the numerical approaches to inverting Eq. (fTfi could be 
adapted to invert Eq. ([2]) to obtain g(r,6). Here we 
choose to adapt the widely used pBasex approaclP and 
expand the desired distribution g{r, 9) as a product of 
basis functions comprising products of Gaussian radial 
functions and Legendre polynomials as angular functions, 


II. INVERSION OF PHOTOFRAGMENT IMAGES WITH 
NON-UNIFORM DETECTION 


"•max ■'max 


*M) = ££ C k ig k i(r, 9), 


k —o i —o 


( 3 ) 


In a photofragment imaging experiment, the Newton 
sphere of charged particles of interest is projected onto 
a 2D position sensitive detector via electrostatic lenses 
which accelerate the charged particles towards the detec¬ 
tor. If the initial distribution is cylindrically symmetric 
about an axis parallel to the detector frame, then the 
3D distribution and its 2D projection are related via the 
Abel transform, 


F(y,z) = 2 



pf(p, z) 
vV -y 2 


dp, 


(i) 


where F{y,z) is the 2D projection, f(p,z) is the 3D dis¬ 
tribution which is assumed cylindrically about the 2 axis, 
and p is the distance to the 2 -axis as illustrated in Fig. |Tj 
Typically, solving the inverse of this equation to recover 
the desired f{p,z) distribution directly is sub-optimal 
due to the sensitivity to experimental noise, and conse¬ 
quently over the years a number of numeric al ap proaches 
have been developed to tackle this problem. 1 S |9|30 | 3i| 
Here we treat the case where the 3D charged particle 
distribution lacks an axis of cylindrical symmetry due to 
a non-uniform detection function. We choose to analyse 
the problem in polar coordinates, which has been shown 
to have advantages in terms of localizing any noise in the 
inversion process to the very centre of the imagep^SUMl 
and also provides a natural description of many optically 
induced processes in atoms and molecules.^ 

We consider the case where we wish to characterize a 
3D cylindrically symmetric distribution g(r,9), which is 
projected onto the detector via a non-uniform detection 
function D(r,9,<j)). Here r, 0, and 4> are respectively 
the radius, polar angle and azimuthal angle describing 


where the basis functions are 

gki{r,6 ) = —^==e 2 J P z ( cos0). (4) 

cry Z7T 

Each radial function has a Gaussian width of a and is 
centred at r k = A r k k where A r k = r max /fcmax and r max 
is the maximum radius of the charged particle cloud con¬ 
sidered. 

The VMI image F(R, 0) can then be written as an 
expansion in the corresponding projected basis functions, 


1c / 
"'max 1 


max ■'max 


F(R, 0) = ££ C k iF k i(R, 0), 


fc =0 i —0 


( 5 ) 


where the projected basis functions are given by 


Fm(R, 0 ) 



rg k i(r, 9)D(r, 0, (f>) ^ 
Vr 2 - R 2 7 ' 


( 6 ) 


In the common case where the projected image is de¬ 
tected on a discrete grid of cartesian pixels, so long as we 
choose the width a of the radial basis functions in Eq. (3) 
to be around 1 pixel, we can express the image Eq. (5) 
in discrete form as 


k 1 

"'max ‘■max 

Fij (Ri, ©,•) = E £ c ^ F ti ( ? ) 

fc=0 1=0 

where i and j are the detector radial and angular pixel 
indices respectively. Here Ri = (i + |)AR is the value 
of R at the centre of pixel (i,j), where AR is the radial 
pixel width, and 0 7 = (j + ^)A0 is the value of 0 at the 
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centre of the pixel where A0 is the angular pixel 

width. The corresponding discretized basis functions are 
given by 


Fii(R i ,e j ) = 


r rg kl {r,0)D(r,9,<t>) 

'Ri \Jr 2 - R 2 


dr, 


where 


cos 9 = 


sin m = 


Ri cos Qj 

r 

Ri sin 0, 


Ri sin 0 j 


r sind r 2 _ r2 cos 2 q. 


( 8 ) 

(9) 

( 10 ) 


The basis functions Eq. (|8j) represent the VMI images 
corresponding to each basis function g k i (r, 6) convoluted 
with the detection function D{r,6,cf)). These projected 
basis functions may be calculated using standard numer¬ 
ical integration methods such as CQUAD in the GSL 
librarj®Tto perform the integration over r in Eq. ([8]). In 
order to obtain the expansion coefficients C k i, the sys¬ 
tem of linear equations represented by the matrix equa¬ 
tion Eq. 0 can then be solved using linear algeb ra tech¬ 
niques such as singular value decompositionPSS We note 
that the method as described is identical to pBasejP in 
the limit D(r, 9, <f >) = 1. 

Once fitted, the C k i coefficients can be used to char¬ 
acterize the charged particle distribution by calculating 
the angular integrated radial spectrum (which is related 
to the speed distribution of the particles) according to 


l[r) 



Y, C k 0 e 2 J 


( 11 ) 


In addition, it is usual to characterize the (radially de¬ 
pendent) angular distribution according to an expansion 
in Legendre polynomials Pi(cos9), 


y/3i(r)Pi(cosO), (12) 

1=0 

where the j3i coefficients are calculated from the C k i co¬ 
efficients as 


f{Q;r) = 


\Z~At: 


Pi(r) 


EiST Cue' 


, r , P2 
r k) 


E fc max _ — 

k =0 


(*-r fc ) 2 


(13) 


alignment and orientation from Coulomb explosion imag¬ 
ing with VMI detection? 13 " 15 17 1E2 

As described in the introduction, in such experiments, 
the molecular sample is first aligned/oriented with strong 
non-resonant laser fields (and sometimes static electric 
fields). Subsequently, in order to measure the degree of 
alignment/orientation produced, an intense probe laser 
pulse with duration much shorter than molecular rotation 
is used to remove multiple electrons from the molecules 
under study. The multiply ionized molecules subse¬ 
quently undergo rapid fragmentation due to Coulomb 
repulsion, and imaging of the resulting ion fragments 
is then used to establish the molecular axis distribu¬ 
tion prior to ionization. Under the assumptions that 
the fragmentation happens rapidly with respect to ro¬ 
tation, and that the fragments recoil in the direction of 
molecular bonds, there is a direct correlation between the 
fragment recoil and the molecular orientation. As men¬ 
tioned previously, the challenge here is to deconvolute 
the non-uniform orientational response of the Coulomb 
explosion process from the measurement in order to yield 
the molecular axis distribution prior to Coulomb explo¬ 
sion. The strategy we propose here is: 

1. Perform a CEI measurement on an isotropic gas 
sample, with the CEI probe polarization direction 
chosen such that an axis of cylindrical symmetry is 
contained in the plane of the detector. 

2. Invert the image from step 1 above to obtain the 
3D distribution of CEI ions by solving Eq. ([ 7 ]) with 
D(r,9,4>) = 1. Since this distribution was obtained 
with an isotropic gas sample, we can obtain the 
orientational dependence of the CEI probe process 
for the CEI probe laser polarization state, pulse 
duration and intensity employed in step 1 above. 

3. Perform a CEI measurement on the 
aligned/oriented molecular sample using the 
same probe polarization state, pulse duration and 
intensity as used in step 1 above. 

4. Invert the VMI image from step 3 using Eq. 0 and 
a detection function derived from step (2) in order 
to deconvolute the orientational dependence of the 
CE process from the observed fragment distribu¬ 
tion, and so obtain the molecular axis distribution. 


which are normalized to /3 0 (r) = 1. 

III. APPLICATION TO MOLECULAR AXIS 
ALIGNMENT AND ORIENTATION PROBED BY 
COULOMB EXPLOSION IMAGING 

A. Methodology 

We turn now to the application of the formalism pre¬ 
sented in Sec. El to the measurement of molecular axis 


We note that we require the molecular axis distribution 
in step 3 above to have an axis of cylindrical symme¬ 
try lying in the plane of the detector in order to apply 
Eq. @. However, step 3 does not require that the same 
geometry of the probe is used as for step 2. For example, 
if using a linearly polarized probe, step 2 requires that 
the linear polarization lies in the plane of the detector, 
but in step 3 we are free to rotate the probe polarization 
to lie in a different direction such as perpendicular to the 
detector plane. For this reason we introduce two frames 
of reference: (i) the detection frame (DF); and (ii) the 
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axis distribution frame (AF). In order to invert the ob¬ 
served image in step 2, we require that the DF possesses 
an axis of cylindrical symmetry lying in the plane of the 
detector in step 1. As such, the detection function will 
have an axis of cylindrical symmetry in the DF. We also 
require that the molecular axis distribution possesses an 
axis of cylindrical symmetry in the plane of the detector 
in step 3. However, in step 3, the DF may be chosen to 
lie in any direction relative to the AF. 

The inversion of the image recorded for the isotropic 
molecular sample in step 2 will yield the fit coefficients for 
the distribution of fragments, C' kl . Here, and in what fol¬ 
lows, we use a prime to denote properties relating to the 
detection function. Under the assumption that the frag¬ 
ments recoil along the direction of the molecular bond, 
this distribution will correspond to the probability of CE 
for each orientation of that bond relative to the laser po¬ 
larization, and so these coefficients can be used to con¬ 
struct the orientational detection function required for 
step 4. 

In the reference frame defined by the detection laser 
polarization, we can write the orientational dependence 
of the CE probe as 

D(9'-r) = cos0'), (14) 

V 47T 


where 9' is measured relative to the cylindrical symmetry 
axis in the DF. The expansion coefficients are cal¬ 
culated from the fit coefficients C' kl according to Eq. (131. 

In order to construct the basis functions Eq. ([8])lor 
step 4, we need to calculate the detection function in 
the AF. The angular dependence of the detection func¬ 
tion expressed in the DF, D(9 r , (j>'\r) is related to the 
angular dependence of the detection function in the AF, 
D(6,(j);r) through a rotation through the Euler angles 
(a, 0 , 7 ) connecting the AF and DF.ES For the present 
case where the detection function has an axis of cylin¬ 
drical symmetry, we can set 7 = 0°. In Appendix [B] we 
show that the detection function in the AF may be ex¬ 
pressed in terms of the (i[{r) coefficients found from the 
probe-alone data inversion as 

D(9, (f>',r) = —7= ft 1 (r)Pi(cos A), (15) 

^ 1 

where 

cos A = cos 12 cos 9 + sin U sin 9 cos(a — </>) 

= cos cos 9 + sin D sin 0(cos a cos <f> + sin a sin cf>). 

(16) 


Eqs. (151 and (161 allow for the evaluation of D(r,9,</)) 
in Eq. (jsj) during the numerical integration over r when 
calculating the basis functions. 

It is important to note that steps 1 and 2 allow us to 
retrieve D(9, </>;r), a detection function dependent upon 
two of the Euler angles ( 9 , <j>) describing molecular orien¬ 
tation in the AF. As such, this detection function is aver¬ 
aged over the third Euler angle x that would be needed to 


specify the molecular orientation.^ This angle describes 
rotation of the molecule around the molecular z-axis. As 
such, this strategy is applicable to extracting the align¬ 
ment/orientation of linear molecules and symmetric ro¬ 
tor molecules. For asymmetric rotor molecules where lo¬ 
calization in x accompanies localization in ^D U7|35| care 
must be taken, and this approach will only apply when 
either the localization in x is small and/or D(9, </, y;;r) 
is independent of %. The latter situation arises for many 
molecules. 


B. Experimental example 

As a demonstration of the approach outlined in Sub¬ 
section ME we apply the strategy to the retrieval 
of the molecular axis distribution of a sample of 1,4- 
diiodobenzene (pDIB) molecules alig ned with a strong 
laser field at 1064 nm wavelength For this molecular 
species, a strong linearly polarized laser field will induce 
alignment of the I-I axis (the most polarizable axis) to¬ 
wards the laser field polarization direction. CEI was used 
to characterize the alignment of the I-I axis, through de¬ 
tection of recoiling I + fragments. Under the assumption 
that the I + fragments recoil axially along the C-I bonds 
in the molecule, the I + recoil direction maps directly to 
the I-I axis direction in the lab frame. For this molecule, 
we expect the dependence of the detection function on 
the angle x to be negligible. 

The molecular sample was prepared in a molecular 
beam with a rotational temperature of ca. 1-2 K. The 
linearly polarized alignment laser field had a pulse dura¬ 
tion of 10 ns which is much longer than the time scale 
for molecular rotation. Consequently, this laser field adi- 
abatically induces molecular axis alignment of the I-I axis 
in the sample, with maxim al alignment occurring at the 
peak of the laser fieldP^ 5 ^ Subsequently a second probe 
laser pulse at 800 nm and with duration of 30 fs was 
timed to arrive at the peak of the alignment laser field. 
This high intensity laser pulse induced Coulomb explo¬ 
sion of the aligned molecules, and the I + ions produced 
were detected with a VMI spectrometer. 

We report here the results of two different studies. In 
the first study a circularly polarized probe laser pulse 
was employed, and the alignment using two different in¬ 
tensities of the aligning laser pulse are compared. In the 
second study a linearly polarized probe laser pulse was 
employed, and we examine the effect of the probe geom¬ 
etry employed by comparing images recorded with the 
probe polarization either parallel or perpendicular to the 
aligning laser polarization. 


1. Circularly polarized probe pulse 

Here we detail experiments carried out with a cir¬ 
cularly polarized laser pulses with intensity of 2 x 
10 14 W/cm 2 . In Fig. [2^, we show the I + VMI data 
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FIG. 2. (a) Experimental I + VMI image recorded for the circularly polarized probe laser alone. The axis of cylindrical 

symmetry (corresponding to the laser propagation direction) lies parallel to the 2 -axis, (b) Corresponding pBasex inverted I + 
image, (c) Radial dependence of the /3i(r) angular parameters corresponding to the pBasex inverted image (solid lines). The 
radial spectrum is also shown (dashed line). 


recorded for the circularly polarized probe alone. Two 
radially separated rings are seen corresponding to two dif¬ 
ferent CE channels. The outermost channel corresponds 
to CE of triply charged pDIB molecules whereas the in¬ 
nermost ring results from CE of doubly charged pDIB 
molecules.® Fig. i> shows the distribution of I + ions 
obtained from the pBasex inversion of the experimental 
VMI image. This inversion was carried out by binning 
the experimental image into a 256x256 polar image and 
solving Eq. Q with fc max = 128, Z max = 10, D(r , 9, <j>) = 1 
(i.e. uniform detection), and a = 1.2 pixels. The coef¬ 
ficients Cu in Eq. (|7l were obtained through projected 
Landweber iteratiorpi® with a projection function set¬ 
ting Cu = 0 if Ck o < 0 at each iteration. Further, due 
to the inversion symmetry of the experiment, only even 
values of l were included in Eq. ( 0 ). 


In Fig. |2}j we show the /3i(r) parameters calculated 
according to Eq. (13), as well as the radial spectrum cal¬ 
culated from Eq. (11). From this plot we can see that in 
regions with non-negligible ion intensity, contributions 
from Pi[r) parameters with l > 4 are negligible, and as 
such the probe detection function is well defined by Cu 


coefficients with l < 4. 


The distribution shown in Fig. 0- represents the de¬ 
tection function for CEI probing with the circularly po¬ 
larized laser pulse for the intensity and pulse duration 
used. It is this distribution that samples the aligned axis 
distribution in the subsequent measurements with laser- 
aligned molecular samples. 

In Fig. 0a) and Fig. 0(b) we show the I + VMI data 
recorded for a sample of molecules aligned with linearly 
polarized laser fields of intensities 1.5 x 10 11 W/crn 2 and 
7.7 x 10 11 W/cm 2 respectively. The aligning laser field 
was polarized parallel to the 2 -axis. This laser field there¬ 
fore induces alignment of the molecular I-I-axis towards 
the 2 -axis. The circularly polarized CE pulse propagated 
parallel to the y-axis such that the light was polarized in 
the X 2 -plane. The observed VMI image therefore corre¬ 
sponds to the molecular axis distribution sampled by the 
detection function of the probe. 


In Fig. 0 and Fig. 01 we show the recovered distribu¬ 
tions of I~ions following deconvolution of the detection 
function determined from the probe alone data follow¬ 
ing the procedure outlined in Subsection III A The in- 
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FIG. 3. (a) and (b): Experimental I + VMI images recorded for laser aligned pDIB probed via Coulomb explosion with a 

circularly polarized laser pulse. The aligning laser polarization is along z, and the probe propagation direction lies along y. 
Images are shown for aligning laser field intensities of (a) 1.5 x 10 11 W/cm 2 , and (b) 7.7 x 10 11 W/cm 2 . (c) and (d): Corre¬ 
sponding pBasex inverted I + images with detection function deconvoluted. Overlaid is a grayscale contour map corresponding 
to the probe alone distribution shown in Fig. In¬ 


versions to recover these distributions were carried out 
with l max = 14 and k max , <r and the number of polar 
bins the same as for the probe-alone data. The detection 
function D(r, 9, <j>) used to construct the basis functions 
(Eq. was calculated from Eq. (141 with A calculated 
from Eq. (16) setting O = 90° and a = 0°. The coeffi¬ 
cients Cm in Eq. |7| were obtained t hroug h Landweber 
iteration with no projection function) 39 * 40 ! Further, due 
to the inversion symmetry of the experiment, only even 
values of l were included in Eq. (J7|. 


Overlaid on each recovered distribution in Fig. n and 
Fig. [3ji is a grayscale contour map corresponding to the 
probe detection function axially integrated over the az¬ 
imuthal angle </>. The calculation of this axially inte¬ 
grated detection function is detailed in Appendix [Cj This 
contour map provides a visual representation of the de¬ 
tection function - its value represents the detection prob¬ 
ability at each value of 9 integrated over all values of </>. 


As can be seen from inspection of Fig. & and Fig. §i, 
the probe detection function samples the molecular axis 
distribution with high efficiency for both aligning laser 
intensities. 

The degree of overlap of the detection function with the 
axis distribution determines the extent to which the full 
molecular axis distribution is sampled, and the reliability 
of the deconvolution process. We can quantify the degree 
of this overlap by evaluating the angular overlap factor 

1 /* 27T n 7T 

0(r) = ---7T / d</> / sin 9 d9f {6; r)D(9, far), 

Anax(r) Jo Jo 

(17) 

where -D max (r) is the maximum value of the angular de¬ 
pendence of the detection function in the AF, D(9 , <f>\ r). 
This integral will take values between 0 (no overlap) and 
1 (maximal overlap). The evaluation of this integral is de¬ 
tailed in Appendix |D| In Fig.[4]we show the radial depen- 
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FIG. 4. Radial dependence of the angular overlap factor corre¬ 
sponding to the pBasex inversion of the images shown in Fig. [3] 
(solid lines) for alignment laser intensities of 1.5 x 10 11 W/cm" 
and 7.7 x 10 11 W/cm 2 , and a circularly polarized probe. The 
corresponding radial spectra are also shown (dashed lines). 


dence of the this overlap factor for the two aligning laser 
intensities employed. The overlap factor is clearly lower 
for the less well aligned distribution at the lower aligning 
laser intensity, reflecting the fact that the broader axis 
distribution extends further into the region of lower prob¬ 
ability of CE by the probe laser pulse, as is also seen by 
comparing Fig. S' and Fig. P- Nonetheless in both cases 
the overlap factor is above 0.9 signifying good sampling 
of the axis distribution. 

In Fig. [5] we show the pi(r) parameters calculated ac¬ 
cording to Eq. dl3l) , as well as the radial spectrum calcu¬ 
lated from Eq. (|11[). For both intensities employed, the 
resulting Pi 4 (r) coefficient remained at 0, and increasing 
the value of ^ max beyond 14 led to no significant change in 
the inverted image. These observations indicate that the 
alignment distribution is well characterized with / max = 
14. At the lower intensity of 1.5 x 10 11 W/cm 2 all Pi co¬ 
efficients are seen to be smaller in magnitude than for the 
higher intensity of 7.7 x 10 11 W/cm 2 , and indeed at the 
lower alignment laser intensity the Pwpr) coefficient was 
seen to be negligible. This is consistent with the high 
alignment laser intensity producing a higher degree of 
molecular axis alignment. In both cases, some large fluc¬ 
tuation in pi (r) values is observed at the largest values of 
r due to the experimental image being slightly truncated 
by the detector edge. 

In Fig.[6]we show the (cos 2 9){r) expectation values for 
the data as well as the radial spectrum calculated from 
Eq. This expectation value is a commonly used 

figure-of-merit for characterizing the degree of molecu¬ 
lar axis alignment. It is important to note that this is 
an expectation value of the molecular axis distribution, 
rather than the commonly used value (cos 2 0), referred 
to as (cos 2 02 d)> which is an expectation value of the pro¬ 
jected image of the axis distribution and which includes 
the effect of the non-uniform detection function. The 


value of (cos 2 0) (r) was calculated acording to 

p 7T ^max 

(cos 2 0)(r) = / ^ Pi(r)Pi (cos 9) cos 2 9 sin Odd. (18) 
■'° ;=o 

We note that for both aligning laser intensities, the 
outermost channel (r ~ 150 pixels) indicates a slightly 
higher degree of molecular axis alignment than the inner¬ 
most channel (r ~ 100 pixels) - this is evidenced by the 
smaller values of the Pi parameters for the inner channel 
compared to the outer channel in Fig. [5j and to a lesser 
extent by the (cos 2 9) expectation values for each chan¬ 
nel. As mentioned, the I + signal in the outermost chan¬ 
nel mainly originates from Coulomb explosion of triply 
charged molecules whereas I + ions in the innermost chan¬ 
nel mainly originate from doubly ionized molecules. The 
triply ionized molecules are produced in the region of 
the probe laser focus where the intensity is highest and 
therefore also in the region where the alignment laser in¬ 
tensity is highest. As such the outermost channel probes 
molecules that are expected to be slightly better aligned 
than molecules probed by the innermost channel. Ad¬ 
ditionally, in the preceding development of our method¬ 
ology we have implicitly assumed that the fragment I + 
ions recoil axially along the direction of the molecular I-I 
axis such that there is a direct correspondence between 
fragment recoil and molecular axis alignment. It is pos¬ 
sible that this axial recoil condition of the I + fragments 
is better fulfilled for the outermost channel than for the 
innermost channel. In general, some deviation from axial 
recoil is expected due to bonding in the multiply charged 
molecular ion created by the probe pulse and to possi¬ 
ble charge-asymmetry in the Coulomb explosion process. 
The effect of the non-axial recoil is expected to slightly 
reduce the degree of measured alignment and it is not 
removed by our deconvolution of the detection function 
- this will be discussed further in a future publicatioiP^. 


2. Linearly polarized probe pulse 

Here we detail a second study carried out with a lin¬ 
early polarized probe pulse of intensity 3.2 x 10 14 W/cm 2 , 
and a linearly polarized aligning pulse of intensity 6.6 x 
10 11 W/cm 2 . In Fig. 0 we show the VMI data recorded 
for the probe laser alone with its polarization along the 
z-axis. The same two CEI channels as observed with 
the circularly polarized probe pulse are evident in the 
VMI data. Fig. 0 shows the corresponding distribu¬ 
tion of I + ions obtained from the pBasex inversion of the 
VMI image. The lower number of counts in this image 
(due to a reduced data collection time) required a more 
coarse binning of the data when carrying out the pBa¬ 
sex inversion compared to the circularly polarized probe 
case. The experimental image was binned into a 128 x 128 
polar image and Eq. 0 was solved with fc max = 64, 
Imax = 10, D(r,0,<p) = 1 (i.e. uniform detection), and 
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FIG. 5. Radial dependence of the /3i(r) angular parameters obtained from the deconvoluted pBasex inversion of the images 
shown in Fig. [ 3 ] (solid lines), (a) alignment laser intensity of 1.5 x 10 11 W/cm 2 . (b) alignment laser intensity of 7.7x 10 11 W/cm 2 . 
The corresponding radial spectra are also shown (dashed lines). 
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FIG. 6. Radial dependence of the (cos 2 8) expecatation val¬ 
ues resulting from inversion of the images shown in Fig. [3] 
(solid lines) for the two alignment laser intensities of 1.5 x 
10 11 W/cm 2 and 7.7 x 10 11 W/cm 2 , and a circularly polar¬ 
ized probe laser. The corresponding radial spectra are also 
shown (dashed lines). 


cr = 1.75 pixels. The coefficients Cm in Eq. 0. were ob¬ 
tained through projected Landweber iteratio n 39 * 40 ^ with 
a projection function setting Cm = 0 if Ck 0 < 0 at each 
iteration. As previously, only even values of l were in¬ 
cluded in Eq. |7|. 


In Fig. [t]i* we show the /3;(r) parameters calculated ac¬ 
cording to Eq. (13), as well as the radial spectrum calcu¬ 
lated from Eq. ©■ As was the case with the circularly 
polarized probe pulse, in regions with non-negligible ion 
intensity, contributions from (r) parameters with l > 4 
are negligible. 


In Fig. E* and Fig. §3 we show VMI data recorded 
for molecules aligned with a linearly polarized laser field 
and probed with a linearly polarized probe pulse in two 
different geometries. For the data in Fig. [8}\ a parallel 
geometry was employed with both the aligning and probe 
laser polarizations along the z-axis. In Fig. i> a perpen¬ 
dicular geometry was used with the aligning laser field 
polarized along the z -axis and the probe laser field po¬ 
larized along the £-axis (perpendicular to the detection 
plane). From these images it is apparent that the relative 
magnitude of the two CEI channels depends strongly on 
the orientation of the probe pulse polarization relative to 
the molecular axis. 


In Fig. [8jc and Fig. [8]l we show the distributions of 
I + ions following deconvolution of the detection function 
determined from the probe alone data following the pro¬ 
cedure outlined in Sec. IIII A1 Overlaid on each recovered 
distribution in Fig. & and Fig. w is a grayscale con¬ 
tour map corresponding to the probe detection function 
axially integrated over the azimuthal angle <f> (Appendix 
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FIG. 7. (a) Experimental I + VMI image recorded for the linearly polarized probe laser alone. The axis of cylindrical symmetry 
(corresponding to the laser propagation direction) lies parallel to the 2 -axis, (b) Corresponding pBasex inverted I + image, 
(c) Radial dependence of the /3z(r) angular parameters corresponding to the pBasex inverted image (solid lines). The radial 
spectrum is also shown (dashed line). 


[0- As can be seen from comparing Fig. |8j: and Fig. [8ji, 
the probe detection function samples the axis distribu¬ 
tion with much lower efficiency when the probe polariza¬ 
tion is along the zr-axis, perpendicular to the direction of 
molecular alignment. Note that in Fig.jsjl, the axially in¬ 
tegrated detection function has a higher efficiency along 
the y-axis than the 2 -axis due to the integration over <f> 
encompassing the zc-axis for 9 = 90°. 

The inversions to recover the distributions in Fig. [8] 
were carried out by binning the experimental image 
into a 256x256 polar image and solving Eq. 0 with 
tax = 128, Z max = 20, and a = 1.2 pixels. The detection 
function D(r, 9 , (f>) used to construct the basis functions 
(Eq. ([8])) was calculated from Eq. (Jl4|) with A calculated 


from Eq. (16) setting 12 = 0° and a = 0° for the paral¬ 


lel geometry and 12 = 90° and a = 0° for the perpen¬ 
dicular geometry. The coefficients Cm in Eq. (pf]) were 
obtained through singular value deconvolution'” which 
was found to give satisfactory results without requiring 
regularization via the projected Landweber iteration. As 
previously, only even values of l were included in Eq. (J7| . 

The radial spectra obtained from the deconvolution 


calculated according to Eq. (11) shown in Fig. [9] clearly 
show that for the perpendicular probe geometry there 
is a reduced relative contribution from the outer CEI 
channel, corresponding to explosion of the triply charged 
parent ion, compared to the doubly charged parent ion 
CEI channel. In addition, other features are observed 
in the perpendicular geometry radial spectra suggesting 
that the relative contributions from different fragmenta¬ 
tion pathways are dependent on molecular orientation.® 
It is interesting to note that these details would not 
be apparent from the raw VMI data before deconvolu¬ 
tion/inversion. 

In Fig. [TO] we show the radial dependence of the an¬ 
gular overlap integral (Eq. ( |T7| ) evaluated as described 
in Appendix [d|) . This plot shows that the perpendicular 
probing geometry has a much lower angular overlap inte¬ 
gral than the parallel probe geometry, and also the circu¬ 
larly polarized probe described in Subsection III B 1 'see 
Fig. [4|. This is due to the fact that for this molecule the 
ionization probability for these CEI channels is greatest 
when the probe laser polarization lies along the I-I molec¬ 
ular axis. However, the ionization probability does not 
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FIG. 8. (a) and (b): Experimental I + VMI images recorded for laser aligned pDIB probed via Coulomb explosion with a 
linearly polarized probe laser pulse. The aligning laser pulse was polarized along the 3 direction. In (a) the probe laser was 
also polarized along the 3-axis. In (b) the probe laser was polarized along the aj-axis (perpendicular to the detector plane), (c) 
and (d): Corresponding pBasex inverted I + images with detection function deconvoluted. Overlaid is a grayscale contour map 
corresponding to the (^-integrated detection function. 


drop to zero when the probe laser is perpendicular to the 
I-I axis. One advantage of the perpendicular probe ge¬ 
ometry is that, for molecules with their I-I axes lying in 
the yz-plane, there is uniform ionization probability with 
respect to molecular rotation about the x axis - as such 
this provides a good measurement of the degree of local¬ 
ization of the I-I molecular axes towards the z-axis for 
those in-plane molecules. It is interesting to note that 
for the perpendicular probe geometry, the inner (dou¬ 
bly charged parent) channel has a higher overlap integral 
than the outer (triply charged parent) channel, showing 
that the inner channel’s dependence on molecular orien¬ 
tation is weaker than that for the outer channel. 


Fig. 11 shows the radial dependence of the (cos 2 9) ex¬ 


pectation value for the parallel and perpendicular probe 
geometries calculated according to Eq. (181. For the 
outer channel the retrieved (cos 2 9) is ca. 0.90 for both 
the parallel and the perpendicular probe geometry. It 


might have been expected that the higher angular over¬ 
lap factor for the parallel probe geometry would lead to 
a higher retrieved value of (cos 2 9) value compared to 
that for the perpendicular probe geometry. On the other 
hand in the parallel geometry the best aligned molecules 
have their I-I axis close to the polarization axis of the 
probe pulse and as such they have an increased proba¬ 
bility of being ionized to higher charge states due to en¬ 
hanced ionizationP^lH'f These higher charged molecular 
ions could fragment into l n+ ions with n > 1 rather than 
into I + , i.e. the best aligned molecules would not lead 
to signal in the I + ion images and would therefore lead 
to a reduced value of (cos 2 9) when determined from the 
I + signal. The almost identical (cos 2 9) value observed 
for the parallel and perpendicular geometries indicates 
that neither the detection overlap factor nor enhanced 
ionization prevents a reliable measurement of the degree 
of alignment for any of the probe geometries using the 
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FIG. 9. Radial spectra obtained from the pBasex inversion 
of the aligned molecule data (Fig. [8| for the parallel (blue) 
and perpendicular (green) probing geometries. The radial 
spectrum obtained from pBasex inversion of the probe-alone 
data for randomly oriented molecules (Fig. [7| is also shown 
(black). All spectra shown are normalized to a maximum 
value of 1. 



FIG. 11. Radial dependence of the (cos 2 6) expectation val¬ 
ues obtained from the pBasex inversion of the data shown in 
Fig. [8] for the parallel (blue) and perpendicular (green) probe 
geometries. The radial spectra are also shown (dashed lines). 


being active for the inner channel. 58 

In Fig. 12 we show the Pi{r) parameters for both the 
parallel and perpendicular probe geometries calculated 
according to Eq. (13). In both cases these coefficients 
show there is negligible contribution from Legendre poly¬ 
nomials beyond 12 for the alignment laser intensity used. 
The Pi{r) coefficients for the two probe geometries agree 
well for the outer CEI channel, but there is a decrease 
in the /3 i(r) coefficients for the perpendicular geometry 
similar to that seen with the (cos 2 9) expectation value. 


IV. CONCLUSION 


FIG. 10. Radial dependence of the angular overlap factor 
corresponding to the pBasex inversion of the images shown in 
Fig. El (solid lines) for the parallel and perpendicular probing 
geometries. The corresponding radial spectra are also shown 
(dashed lines). 


algorithm presented here. In addition, it is clear that 
although the relative weightings of different fragmenta¬ 
tion channels depends on the molecular orientation (as 
seen from the radial spectra, Fig. [9])P1 this is correctly 
accounted for in the retrival algorithm presented, as evi¬ 
denced by the consistent (cos 2 6) values for the two probe 
geometries. For the inner channel, the retrieved (cos 2 9) 
value is lower for the perpendicular probe geometry com¬ 
pared to the parallel probe geometry. Since both probe 
geometries are sampling an identical molecular axis dis¬ 
tribution, the retrieved (cos 2 9) value should be the same 
in both cases, as is observed for the outer channel. As 
with the slightly reduced degree of alignment observed 
for the inner channel when probing with the circularly 
polarized probe in Subsection |IIIB 1[ we attribute this 
difference as arising due to non-axial recoil geometries 


We have proposed a method for deconvoluting a non- 
uniform detection function from velocity map imaging 
experiments provided the detection function is measur¬ 
able independently. Experimentally we demonstrated 
this technique by recovering the axis distribution of ID 
aligned pDIB molecules using laser-induced Coulomb ex¬ 
plosion imaging. A major advantage of the technique 
is that it allows a transferable and complete character¬ 
ization of the axis distribution of aligned molecules. In 
particular, (cos 2 9) can be determined. This represents a 
measure of the true degree of alignment rather than the 
usual (cos 2 02 d) value, determined directly from 2D ion 
images, which is strongly biased by the orientational de¬ 
pendence of the probe process. Furthermore, the method 
also provides higher moments of the axis distribution 
and as such a complete characterization alignment of the 
molecules is possible. 
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FIG. 12. Radial dependence of the /3i(r) angular parameters obtained from the deconvoluted pBasex inversion of the images 
shown in Fig. [8] (solid lines), (a) parallel probe polarization geometry (b) perpendicular probe polarization geometry. The 
corresponding radial spectra are also shown (dashed lines). 


Appendix A: Abel inversion in spherical polar coordinates 


Noting that 


For a specific value of r (proportional to the particle 
velocity), the relationship between the 3D distribution of 
detected particles, /(r, 0,</>) = g(r,9)D(r,9,(j)) : and the 
observed 2D projected image F(R,Q) is 


F(R, 0; r)S Re = f(r , 0, (j>)S e(t> 


(Al) 


where S R q is the elementary surface on the detector, and 
Sec/, is the elementary surface on the sphere of radius r, 


(A2) 

(A3) 


S R e = R di?.d0, 

Se<p = r 2 sin 9 d 9 d<j). 

Srq may be re-written as 

= A|J|d0d0, (A4) 

where the determinant of the Jacobian J is given by 

8R 50 _ 90 OR 
89 d(p 89 8(j) 


IJI = 


(A5) 


Subs titu ting Eq. (A3), Eq. (A4) and Eq. (A5) into 
Eq. dAlj) we can write 


F(R, 0; r) = 


f(r, 9)r 2 sin0 

W\ ■ 


(A6) 


R = r\J cos 2 9 + sin 2 9 sin 2 <f>, 
0 = arctan(sin (p tan 9 ), 


we can evaluate the Jacobian in Eq. (A5) as 
2 r sin 2 (0) cos{4>) 


Ul = 


a/2cos(2 9) cos 2 (i i>) — cos(2 4 >) + 3 
Using the relationships 

'tan© 


= arcsm 


0 = arccos 


tan(? 
r cos 9 
R 


we can subsitute Eq. (A8) into Eq. (A 61 to obtain 


F(R, 0; r) = 


rf{r,6,(j>) 
'Jr 2 — R 2 


with (f> given by Eq. (A91 and 9 given by 

. R cos 0 s 

9 = arccos 


(A7a) 

(A7b) 


(AS) 

(A9) 

(A10) 

L 

(All) 

(A12) 
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Since in general we have more than a single kinetic energy 
present in the 3D distribution, we have to integrate over 
r > R in order to calculate the projection intensity at 


F(R,S) 



r fir, 0,4) 

Vr 2 - R 2 


(A13) 


Appendix B: Rotation of the the detection function from 
the DF to the AF 


The detection function in the DF (Eq. @) may be 
re-written as an expansion in spherical harmonics as 


D(f,V;r) = 7 C^ 


i=0 


47T 

21 + 1 


Pl{r)Y l0 {6', <!>')■ (Bl) 


The detection function in the AF is related to the detec¬ 
tion function in the DF through rotation by the Euler 
angles (a, fl, 7 ).^ The detection function in the AF can 
be written as 


D{0,4\ r) 



E 


1=0 


47T 

21 + 1 


m 


m=l 

X E 

m——l 


(B2) 


where 1, 7 ) are the Wigner rotation matrices. 

Expressing the rotation matrix D l m0 (a,fl 1 7 ) in terms of 
a spherical harmonic yields 


D(0,4;r) 



E 


1 -0 


47T 

21 + 1 




m=l 

X 53 Y; m (n,a)Y lm (0,4), 

m——l 


(B3) 


and 


D' 0 (a,D, 7 )=P i ( cosF!), (C2) 

we can evaluate the axially integrated detection function 
as 


n2n 

D(0\ r) = / D(0,(j)-,r)d<j) 

Jo 


= V^E P 1 (r)Pi (cos Q)Pi (cos 9). 


(C3) 


1=0 


Appendix D: Overlap function evaluation 


The angular distribution of molecular axes in Eq. (12) 
can be re-expressed in terms of spherical harmonics as 

= < D1 > 

Substituting Eq. (B2| and Eq. (Dll into Eq. gives 


-1 •'max / -1 

° (r) = ^W?)gV«TI' 3 ' (r) 




X ^2\l ^F+l 13 ' 1 '^ E D rn' o( a ’ Q ’7) 


i'=l' 


(D2) 




'=0 
/*27T /*7T 

x / / Yio(9,4)Yi'm'{0,4) sind d 0 dip, 

Jo Jo 

Evaluating the integral us ing t he orthogonality of spher¬ 
ical harmonicJ^ and Eq. (C2) gives 


0(r) = D 1 ^ E y^Pi( r )Pl( r ) p i( C0BS1 )- (D3) 


The product of two spherical harmonics can be con¬ 
tracted by the spherical harmonic addition theoremp 4 ! 


4-7T 

2Z + 1 


5] Y* m [Q,, a)Yi m (9,4) = Pi {cos A), (B4) 


1 = — l 


where A is g iven by Eq. (16 1 . Substitution of Eq. (B4) 
into Eq. (B3) gives Eq. (15) 


Appendix C: Detection function integrated over (j> 


In order to visualize how the detection function sam¬ 
ples the axis distribution, it is h elpf ul to calculate the 
detection function in the AF Eq. (B2) integrated over cj>. 
Noting that 



Yim{9,4>) # = 27t 


2Z + 1 
47T 


Pi(cos9)S m0 , 


(Cl) 
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